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I. INTRODUCTION 



A major goal of nuclear structure theory is to make quantitative calculations of low-energy 
nuclear observables starting from microscopic inter-nucleon forces. Renormalization group 
(RG) methods can be used to soften the short-range repulsion and tensor components of 
available initial interactions so that convergence of nuclear structure calculations is greatly 
accelerated jllH]. A major complication is that these transformations change the short-range 
many-body forces. In fact, any softening transformation will induce many-body interactions 
in the course of renormalizing the matrix elements in a lower sector. To account for these 
changes, we must include consistently evolved three-body (and possibly higher) forces in 
structure calculations. 

A previous letter [3| presented the first such evolution of three-body forces in free space 
by using the Similarity Renormalization Group (SRG) 0-01 • The SRG offers an approach 
to evolving many-body forces that is technically simpler than other unitary RG formula- 
tions. Irrespective of the chosen initial Hamiltonian, the evolution produces a variational 
Hamiltonian and enables smooth extrapolation of results, in contrast to Lee-Suzuki 10] type 
transformations which produce results that are model-space dependent (in both iV max and 
A) [llj]. While the SRG induces many-body forces as a product of renormalization, these 
terms come in a hierarchy of decreasing strength if a hierarchy is initially present. Partic- 
ularly useful in an analysis of such a hierarchy are chiral effective field theories (xEFTs), 
which provide a systematic construction of many-body forces as the initial input to our 
evolution calculations. [12]. Our results expand on prior evidence that the SRG explicitly 
preserves the initial EFT many-body hierarchy as it improves convergence properties of 
evolved Hamiltonians. 

Section HT1 reviews some background material on how the SRG is applied in these calcula- 
tions. In Section UTTl we explore the convergence properties of the renormalized Hamiltonians, 
including new A=6 calculations. In Section IIVI we present the calculations as a function of 
the evolution parameter, and explore the effect of SRG flow on other initial interactions. 
Section |V] dives deeper into the analysis of how the SRG acts to evolve the input interac- 



tion, expanding upon the analysis done for one-dimensional models [13J. We make a brief 
advertisement of operator evolution and conclude with comments on the future use of this 
approach. 



II. BACKGROUND 

As implemented in Refs 0, [1] for nuclear physics, the SRG is a series of unitary transfor- 
mations, U\, of the free-space Hamiltonian, 

H x = U X H X=00 U{ , (1) 

labeled by a momentum parameter A that runs from oo toward zero, which keeps track of 
the sequence of Hamiltonians (s = 1/A 4 is also used elsewhere 0,11]). These transformations 
are implemented as a flow equation in A (in units where h 2 = M = 1), 

^ = 4[P\*],*d, (2) 

whose form guarantees that the Hx's are unitarily equivalent 0, 0]. Once the Hamiltonian 
has been evolved we also have the option to build the unitary transformation operator 



2 



directly as a sum over outer products of the evolved and unevolved wavefunctions: 



i/a = X>«(a)W«(o)|, (3) 

a 

where a is an index over the states in the chosen configuration space. This feature is useful 
in applications to external operators [HI]. Note that U\ can also be evolved directly and the 
choice of method is open to efficiency and convenience for a particular use. 

The appearance of the nucleon kinetic energy T in Eq. ([2]) leads to high- and low- 
momentum parts of H\ being decoupled, which means softer and more convergent po- 



tentials [15j . This is evident in a partial- wave momentum basis, where matrix elements 
(k\H\\k') connecting states with kinetic energies differing by more than A 2 are suppressed 
by e - ^' ~ h ' ) / A factors and therefore the states decouple as A decreases. However, decou- 
pling also results from replacing T in Eq. fl2]) with other operators @, 0, [l6, 17]. The optimal 



range for A is not yet established and also depends on the system. Previous experience with 
SRG and other low-momentum potentials suggested that running to about A = 2.0fm _1 is a 
good compromise between improved convergence from decoupling and the growth of induced 
many-body interactions [l5|. Current results show that this limit might be extended as far 
as A = l.Ofm -1 , at least for lighter nuclei. 

One formal way to see how the two-, three-, and higher-body potentials evolve is to 
decompose H\ in second-quantized form 0. We can write a general A-body Hamiltonian 
as: 

Hx = J2 T V a i a 3 + 2j2 V W,A a ) aiak + ^2 Yl V mrnn,X a WA a narnai + • • • , (4) 

ij ijkl ijklmn 

where a\ and at are creation and destruction operators with respect to the vacuum in some 
single-particle momentum basis. The quantities T^, V^, and vM lmn denote matrix ele- 

ments of their respective operators. Equation (j4j) defines Ty, V^y A , ^W JmnA , ...as the 
one-body, two-body, three-body, . . . matrix elements at each A. By evaluating the commu- 
tators in Eq. (j2J) using H\ from Eq. (j3J), and normal ordering the resulting terms of cre- 
ation/annihilation operators, we find that higher-body potentials are generated with each 
step in A, even if initially there are only two-body potentials. We note that in this paper we 
are not actually evolving in a single-particle basis as indicated in Eq. (j3J) , but nothing a pri- 
ori prevents it as a choice of basis. In particular, the center-of-mass solutions will factor out 
in the properly truncated Harmonic Oscillator single-particle basis. Furthermore the SRG 
will not mix different center-of-mass solutions since Tqm commutes with the Hamiltonian. 

Here we are normal ordering with respect to the vacuum, as opposed to the in-medium 
SRG which normal orders with respect to a non-vacuum reference state. With in-medium 
normal ordering, SRG evolution generates an A-dependent rearrangement of the higher- 
body contributions to the evolved interaction; the density-dependent 0-, 1-, and 2-body 
normal-ordered interactions are found to absorb the dominant free-space many-body inter- 
actions [HI . For free-space normal ordering, matrix elements in a given sector are determined 
completely by evolution in that sector. In addition, each A-body sector contains as a subset 
the (A — l)-body sector evolutions. Thus, when applied in an A-body subspace, the SRG 

(2) 

will "induce" A-body forces, with (T) fixed, (V^ ) determined completely in the A = 2 

rn\ (n\ (2) (-1) 

subspace with no dependence on (V" A ), (V^ ) determined in A = 3 given (V^ ) and (V^ =0 ), 
and so on. 
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Because only the Hamiltonian enters the SRG evolution equations, there are no difficulties 
from having to solve T matrices (of the Lippmann-Schwinger equation) in all channels for 
different A-body systems [If. However, in a momentum basis the presence of spectator 
nucleons requires solving separate equations for each set of (V^ ) matrix elements. In 



Ref. [19], a diagrammatic approach was introduced to handle this decomposition. But while 
it is natural to solve Eq. ([2]) in momentum representation, it is an operator equation allowing 
us to use any convenient basis. Here we evolve in a discrete basis of Jacobi-coordinate 
harmonic oscillator wave functions, where spectators are handled without a decomposition 
and induced many-body interaction matrix elements can be directly identified. Having 
chosen a basis, we obtain coupled first-order differential equations for the matrix elements 
of the flowing Hamiltonian H\, where the right side of Eq. ([2]) is evaluated using simple 
matrix multiplications. 

The procedures used here build directly on Ref. [HI], which presents a one- dimensional 
implementation of our approach along with a general analysis of the evolving many-body 
hierarchy. We start by evolving H\ in the A = 2 subsystem, which completely fixes the two- 

(2) 

body matrix elements (V A ). Next, by evolving H\ in the A = 3 subsystem we determine 
the combined two-plus-three-body matrix elements. We can isolate the three-body matrix 
elements by subtracting the evolved (V^) elements in the A = 3 basis [l3|. Having obtained 
the separate NN and NNN matrix elements, we can apply them unchanged to any nucleus. 
We are also free to include any initial three-nucleon force in the initial Hamiltonian without 
changing the procedure. If applied to A > 4, four-body (and higher) forces will not be 
included and so the transformations will be only approximately unitary. The questions to 
be addressed are whether the decreasing hierarchy of many-body forces is maintained and 
whether the induced four-body contribution is unnaturally large. We summarize in Table |T] 
the different calculations to be made here for 3 H, 4 He, and 6 Li to confront these questions. 
These calculations will also be made when other nuclei are considered. 



TABLE I: Definitions of the various calculations. 
NN-only — No initial NNN interaction 

and do not keep NNN-induced interaction. 
NN + NNN-induced — No initial NNN interaction 

but keep the SRG-induced NNN interaction. 
NN + NNN — Include an initial NNN interaction 

and keep the SRG-induced NNN interaction. 



Hamiltonians obtained via free-space SRG evolution are independent of the basis choice. 
Up to truncations induced by conversion to a particular basis, a Hamiltonian evolved to a 
given A reproduces the results of a Hamiltonian evolved to the same A in a different basis. 
Two types of truncations, in model-space size and A, are relevant to controlling the quality 
and consistency of SRG evolved interactions. Our calculations are performed in the Jacobi 
coordinate harmonic oscillator (HO) basis of the No-Core Shell Model (NCSM) [2(j. This is a 
translationally invariant, anti-symmetric basis for each A-body sector, in which a complete 
set of states in the model space defines the maximum excitation of N max h~Q above the 
minimum energy configuration, where Q is the harmonic oscillator parameter. Hamiltonians 
are derived and evolved in this basis and then switched to a Slater determinant basis. The 
Jacobi coordinates used to build this basis have a convenient normalization that treats 
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all A-body clusters on an equal footing. Operators in an A-body space, like the A-body 
Hamiltonian, can be embedded in an (A + n)-body space in a straightforward manner. Due 
to the antisymmetric nature of the basis, they need only be multiplied by a combinatoric 
factor: 

'A + n\ _ (A + n)\ 

A J ~ A\(n)\ ' [ ' 

For example, a 3-body system has Q) = 3 pairs, a 4-body system has ( 2 ) = 6 pairs and 
(3) = 4 triplets, etc. This embedding factor was a direct predictor of behavior in one 
dimension [13j, but in the realistic case, many physical constraints may complicate the end 
results including, but not limited to, Pauli blocking, angular momentum selection rules, and 
cancellations intrinsic to the initial Hamiltonian. 

A major drawback of the HO basis is its single intrinsic scale, Ml, which is problematic for 
systems with multiple scales. However it is a widely used basis in part because it facilitates 
the separation of spurious center-of-mass solutions and vital to the translationally invariant 
physics of nuclear structure calculations. To understand the cutoffs inherent in the finite 
oscillator basis we can consider that the local maxima in the harmonic oscillator function 
are essentially Gaussians modulated by polynomial terms up to iV max . These maxima, in 
the momentum and coordinate space representations, will be correlated with the high and 



low-momentum cutoffs respectively. These cutoffs have the large iV max behavior: [21 



A uv ~ a/ mN max htt and A m ~ J— — . (6) 

1 \ xr 



mMl 



" max 



When Ml grows large, individual oscillations are large and lose resolution on the small details 
in the momentum basis potential that correspond to large r structures. However, high N max 
polynomials have many small oscillations at low momenta compensating for the large Ml 
value. Thus, Air is lowered and Auv is raised by increasing A" max as expected when the 
basis is extended towards completeness. Note that only the value of N mscK , not Ml, affects 
the number of matrix elements in the basis, so the computational cost is the same for each 
Ml. Changing Ml effects the balance between Air and Auv completeness. 

Note that the behavior attributed to Air does not manifest as an explicit cutoff in the 
momentum representation, but rather a distortion of matrix elements at low momentum. 
Specifically, the effective cutoff operator in momentum representation displays bands of 
ringing artifacts along the off-diagonal direction that ultimately behave as a cutoff; both 
smaller Ml and larger A" max bases alleviate this effect, as is apparent from Eq. (jH]). 

Because of computational constraints we were forced to apply separate truncations, 
AA2max and A^A3maxj to the A = 2 and A = 3 sectors of the initial Hamiltonian (see Table HT1) . 
In previous work [3j with 3 H and 4 He and an initial xEFT interaction, we found that A^max 
= A"A3max = 32 was sufficient because these nuclei are less sensitive to the asymptotic behav- 
ior of the oscillator wavefunctions. However, for 6 Li or 4 He using a harder potential (such 
as Argonne Vis), larger space was required for the initial NN Hamiltonian. When needed 
we used A^max = 300, which is more than enough to accommodate any potential at any 
relevant Ml. In some calculations we were restricted to A^max = 196, but these cases are 
also converged to the keV level. The A = 3 basis size grows much more quickly so that 
the evolution of Hamiltonians above A^max — 40 are very intense computations. In our 
final results there is a slight effective truncation of the induced three-body forces, but this 
is only a truncation of the initial interaction that is then evolved. Additionally, the NNN 
interaction is a perturbative correction to the NN so this truncation has a small impact on 
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our results. These truncation issues are addressed throughout the paper when discussing 
the convergence properties of our results. 



TABLE II: Definitions of the truncations used on initial Hamiltonians. 
N m3X : The size of the final A-body space 
AA3max : The size of the basis for initial three-body matrix elements 
AA2max : The size of the basis for initial two-body matrix elements 



The present calculations make use of both the Jacobi oscillator basis described above 
and a Slater determinant oscillator basis often referred to as the m-scheme. The size of 
the Jacobi basis scales well with N max but poorly with A due to the effort involved in 
antisymmetrization. A Slater determinant basis trades ease of antisymmetrization for very 
large (dimensions into the billions) but sparse matrices, solvable with the Lanczos algorithm. 
Given the convergence advantages provided by the SRG, we obtain the initial Hamiltonians 
in the Jacobi basis, evolve them in A=2 and 3, and transform them into a Slater determinant 
basis for use in existing configuration interaction (CI) codes. 

Our calculations are limited by the size of the input three-body interaction file: the 
present code is not able to distribute the matrix elements among nodes and therefore must 
hold the entire A-body Hamiltonian on each node. For our calculations of 6 Li at A" max = 8, 
the 3-body file is 13Gb of matrix elements in addition to 2-body matrix elements already 
stored. This is the largest calculation possible on most available nodes with 16 or 32 Gb 
of memory; at N max = 10 the 3-body file is 33Gb. However, the m-scheme code MFDn of 
Vary et al. j22[ is capable of distributing the input matrix elements among several nodes and 
efforts are underway to perform these calculations in larger spaces. Furthermore, the size of 
these input files can be dramatically reduced in the future by implementing a compression 
scheme based on angular momentum couplings and calculations could be extended with 



the importance truncation method [23] , |24j. Future calculations of 8 Be, 10 B, and 12 C are 
planned. Recent results of 12 C and 16 using importance truncation provide a benchmark 
for future efforts [25| . 

For the lower A's in the lighter nuclei, our predictions for ground-state energies are fully 
converged. However, in other cases we need to extrapolate the energies to A" max = oo. Here 
we use the same extrapolation procedure applied in [2J. The model used for ground-state 
energies is 

E a i = -Eoo + A a e , (7) 

where a labels the Ml values, i the Af max values for each a, and A a and h a are constants. 
The goal of a fit to the following calculations is to determine the common parameter E^, 
which is the estimate for the ground-state energy extrapolated to A" max = 00. 

This can be cast as a one-dimensional constrained minimization problem with the function 

3 (£00) = ^(log(£ ai - £00) - a a - b a Ni) 2 /a 2 ai , (8) 

where the {a a } and {b a } are determined directly within the function g by invoking 
a constrained linear least-squares minimization routine. The constraint is the bound 
-Eqo ^ min({E ai }), where E^ < and "min" means "most negative". We can also allow 
for weights depending on A" max and/or Ml In the present investigation, we apply Eq. ([8]) 
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with only the hVt value that yields the lowest energy in the largest space, weighting different 
-/Vmax by the slope of the energy vs. iV max , and using the spread of results from neighboring 
Ml values to determine a conservative confidence interval for the extrapolation. Alternative 
approaches to extrapolation in the NCSM are described in Ref. [26| . 

We have considered a variety of interactions as initial inputs to the SRG evolution, 
including chiral EFT, Argonne Vi§ 27J, and CD-Bon n |28l| . The initial (A = oo) chiral NN 
potential is the 500 MeV N 3 LO interaction from Ref. [29]. With the chiral potential we also 
have available an initial NNN potential at N 2 LO in the local form of Ref. j3l[, with 
constants cd = —0.2 and ce = —0.205 fit to the average of triton and 3 He binding energies 



and to triton beta decay as described in Ref. [32J. NCSM calculations with these initial chiral 
interactions and the parameter set in Table I of Ref. j32[ yield energies of —8.473(4) MeV 
for 3 H and -28.50(2) MeV for 4 He compared with -8.482 MeV and -28.296 MeV from 
experiment, respectively. There is a 20keV uncertainty in the calculation of 4 He from 
incomplete convergence but a 200 keV discrepancy with experiment. The latter is consistent 

These provide 
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with the omission of three- and four-body chiral interactions at N 3 LO 
the scale for assessing whether induced four-body contributions are important compared to 
other uncertainties. The best result for 4 He using the AV18 potential is —24.23(1) MeV [HJ 
and using CD-Bonn we compare to —26.1(1) MeV. Here there are larger discrepancies with 
the experimental values due to the omission of consistent initial three-body interactions, but 
these calculations are still useful to assess the effects of induced NNN. 

For 6 Li calculations we use only the chiral interactions at N 3 LO [29[ for NN and at N 2 LO 
for NNN in the form described above. The best existing binding energy with the N 3 LO 
interaction, using a Lee-Suzuki based renormalization up to iV max = 14, is 28.5 ± 0.5 MeV. 
With NNN included the converged value is 32.5±0.5 MeV [35]. The truncations analogous to 
-VA2max and A^A3max for these calculations were A^A2max = 400 and A^A3max = 40, equivalent to 
the initial Hamiltonian inputs for the present work. Results are generally not dependent on 
the particular values of the LECs in a range of cd ~ —2 to + 2 (with cd and ce constrained 
by the fit to 3 H 32] ) but some observables may be particularly sensitive as discussed in 
Ref. S. 



III. CONVERGENCE 

In Fig. [H we show the triton ground-state energy as a function of the oscillator basis 
size, iV max . The convergence of the bare interaction is compared with the SRG evolved to A 
= 2.0 and 1.5 fm _1 . The oscillator parameter HVL in each case was chosen to optimize the 
convergence of each Hamiltonian. We also compare to a Lee-Suzuki (LS) calculation (green 



squares), which has been used in the NCSM to greatly improve convergence [36|, l37j. All of 
these effective interactions result from unitary transformations. The LS is done within the 
model space of a target nucleus, in contrast to the free-space transformation of the SRG, 
which yields nucleus-independent matrix elements. Consequently, the LS results are non- 
variational independent calculations at each N m3iX while the SRG-evolved Hamiltonians can 
be simply truncated to produce the curves shown. A dramatic improvement in convergence 
rate compared to the initial interaction is seen even though the xEFT initial interaction is 
relatively soft. The SRG acts to decouple high-momentum degrees of freedom so the UV 
part converges faster with respect to iV max . Thus, once evolved, a much smaller N max basis 
is adequate for a particular accuracy. 

Figure [2] illustrates for 4 He the same rapid convergence with A^ax of an SRG-evolved 
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FIG. 1: (Color online) Ground-state energy of 3 H as a function of the basis size iV max for an N 3 LO 
NN interaction [2^] with an initial NNN interaction 12, 32]. Unevolved ("bare") and Lee-Suzuki 
(L-S) results with Ml = 28 MeV are compared with SRG evolved to A = 2.0 fm -1 and 1.5£m — 1 
with Ml = 20 MeV. 



interaction. However, in this case the asymptotic value of the energy differs slightly because 
of the omitted induced four-body contribution. The difference can be as large as 100 keV 
for A = 1.0 but no larger than 50 keV for the substantially evolved A = 2.0 fm - . The SRG- 
evolved asymptotic values for different Ml (solid vs. dotted curves) differ by only lOkeV, 
so the gap between the converged bare/L-S result and the SRG result is dominated by the 
induced NNNN rather than incomplete convergence. Convergence is even faster for lower A 
values, ensuring a useful range for the analysis of few-body systems. However, because of 
the strong density dependence of four-nucleon forces, it will be important to monitor the size 
of the induced four-body contributions for heavier nuclei and nuclear matter. In Section [V] 
we present a tool for analyzing the growth of induced many-body forces. 

Also evident in Fig. [2] is the evolving dependence on Ml. Calculations are variational 
in Ml with the optimal value indicating balance between Kuv and h-m- For the initial 
Hamiltonian, the limit of iV max =20 is small and the larger Ml is necessary to provide a 
sufficient Ayy. The larger IR cutoff due to higher Ml is less of a problem than the smaller 
UV cutoff due to the low iV max . However, if the initial Hamiltonian is evolved in the iV max =32 
space, then more UV information is shifted down into the A r max =20 space. Now the high IR 
cutoff is more significant and a lower Ml is more optimal. In the figure, one can see that the 
Ml = 28 calculation (solid curves) has significantly better convergence properties for lower 
A, especially at the small N max (< 8) that is crucial for larger A. 

Evolving Hamiltonians such as CD-Bonn and Argonne Vis also results in much improved 
convergence properties, as seen in Fig. [31 Here a large initial A=2 cutoff, A?A2niax> is cru- 
cial due to the strong high-momentum components in the AV18 potential. However, at 
-^A2max=196 the NN-only results are converged with respect to variation in A^max to within 
IkeV. For CD-Bonn we found that N^ max —40 and iV max =20 was sufficient to converge re- 
sults to within 30 keV. Calculations for the Argonne potential require a bit more effort, 
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FIG. 2: (Color online) Ground-state energy of 4 He as a function of the basis size N max for an N 3 LO 



NN interaction 29f] with an initial NNN interaction 
with Lee-Suzuki (L-S) and SRG evolved to A 
black arrow indicates the experimental value. 



12 



32] 



Unevolved (A) results are compared 
2.0 and 1.2 fin -1 at Ml = 28 and 36 MeV. The 



© — @ y _ 


1.0 fm" 1 




@ — n ^ 


1.2 fm" 1 






1.5 fm"' 






2.0 fm"' 






3.0 fm"' 






N =40 

A3max 



12 14 16 

N cut 



> 



M 

C 
W 

i 

-a 

C 
P 
O 
Sh 

a 



-24.5 



-25.5 



-26.5 



© © T| 


1.0 fm" 1 


O H 


1.2 fm" 1 




1.5 fm" 1 




2.0 fm" 1 




3.0 fm" 1 




10 



12 14 16 
N cut 

max 



18 20 22 



FIG. 3: (Color online) Ground-state energy of 4 He for select A as a function of the basis size N n 
for the AV18 [27( and CD-Bonn [28| interactions. Results are shown for Ml = 44 MeV with iV A 3 n 
= 40. Dotted lines indicate current best results for these potentials 34]. 



obtaining convergence with respect to A^max to within 130 keV. The optimal frequency for 
both of these interactions evolved to A = 2.0 was found to be Ml = 44 MeV. Larger values 
for A^A3max are possible for these potentials, but the current level of convergence is sufficient 
to observe the qualitative behavior of the SRG in Section [TVl 

A multitude of chiral EFT interactions are available for use in initial Hamiltonians. We 
have used here the one version (500 MeV N 3 LO from Ref. 29J) for which the accompanying 
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three-nucleon terms have been rigorously fit to data. An important task in the future will 
be to apply SRG techniques to many more available interactions to compare and contrast 
them. Some of these can be significantly harder than that chosen here, so running to low A 
would be especially important though computationally expensive. 




T 



T 



T 



— \= 1(X) I'm 
M X = 5.0 fin" 
* * X = 3.0 fin" 
" X = 2.0 fm" 
U = 1.5 fm" 
^X = 1.0 fm" 

NN+NNN 

Ml = 20 MeV 



N . . =300 

Azmax 

N =40 

A3max 



IT 



1L 



X 



10 12 
N cut 

max 



14 16 



18 



FIG. 4: (Color online) Ground-state energy of 6 Li as a function of the basis size iV max for an initial 



N 3 LO NN interaction [29( and an initial N 2 LO NNN interaction [12|,|32|] evolved to various A. Note 
the large scale and complete lack of convergence for unsoftened Hamiltonians (A = 100 fm -1 ). The 
dotted line is the best Lee-Suzuki result. 



Figure H] shows just how important a softening transformation is to achieve convergence in 
light nuclei. Here we plot 6 Li binding energies up to A" max = 8 for several A's from 100 to 1.0 
fm" 1 . A meaningful extrapolation is simply not possible, even with the relatively soft chiral 
potential, without some form of softening renormalization like the SRG or Lee-Suzuki type 
transformations. A key advantage of the SRG program is the ability to perform systematic 
extrapolations to spaces that are computationally inaccessible. 

Compare this to the case of 4 He where the initial chiral EFT Hamiltonian is sufficiently 
soft to produce a nearly converged result at A" max = 20. For A = 6 we are restricted to 
smaller N m & x because the basis scales with N max much faster than at A = 4. In addition the 
radius of 6 Li is larger and requires a lower IR cutoff. Thus more oscillator basis states of 
the initial interaction are required to accurately describe this nucleus. In other words, even 
if we could perform the N max = 20 calculation of 6 Li, it would still not converge as well as 
4 He does at that level. With the SRG the information of the larger basis can be moved into 
a smaller space in a smooth controlled way. 

In Fig. [5] we show the convergence patterns of the ground state of 6 Li on a more detailed 
scale. Here the three different calculations of Tabled] are shown side-by-side for clarity, with 
several A values. The initial Hamiltonian was defined by the truncations A^max — 300 and 
A^ASmax — 40 and the truncation errors from (or convergence with respect to) these limits 
are 1 keV and 80 keV respectively. The 6 Li calculation was performed up to A r max =8 for 
the three-body matrix element versions, and A^ max =10 for NN-only (possible with only two- 
body matrix elements). Further calculations with three-body matrix elements will require 
a distributed memory approach like MFDn of Vary et al [22] and the other techniques 
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FIG. 5: (Color online) Ground-state energy of 6 Li as a function of the basis size iV max for an 
initial N 3 LO NN interaction [29] and an initial N 2 LO NNN interaction 12, 32] evolved to various 
A. Here the initial NN potential was included up to iVA2max = 300 and the NNN up to N^max = 
40. Results are compared with the best Lee-Suzuki shown by the thin black dotted line (see text). 
The dotted lines extending from the curves are examples of the extrapolations made throughout. 



mentioned above; such codes will be used for future calculations of larger p-shell nuclei. 
The straight dotted line shows the converged Lee-Suzuki results for NN-only and NN+NNN 
calculations. 

While we can see the improvement in convergence by the data points alone, we can only 
measure the effect of induced many-body forces by considering the extrapolated converged 
value for each A. We provide a sample of the extrapolations we will use to assess the 
converged values for 6 Li, in the form of the dotted lines extending from each curve. The 
spread of these lines at large iV max is the chief indicator of remaining scale dependence in 
the results. In the NN-only plot on the left, one can clearly see the large spread due to 
omitted induced three-body matrix elements. This spread is decreased in the other plots 
by including 3-body matrix elements, first in the center by including those induced by the 
renormalization, and on the right by also including the initial NNN strength. In both NNN- 
inclusive plots the curves have a smooth qualitative progression from higher to lower A, 
indicating less interplay between attractive and repulsive components of the interaction, as 
discussed in Section IVl 

The size of induced 4- to 6-body forces in this calculation is estimated by measuring the 
spread of the lines in the center and right plots, or alternatively by considering the slope of 
the binding energy as a function of lambda. However, the spread is actually smaller than it 
appears here because only the A=1.5 and 1.8 curves are satisfactorily converged at this Ml. 
The curves for other A's are optimal at different values of Ml and their converged values are 
not accurately represented in this simple example. However, the full extrapolation procedure 
does indeed take this Ml dependence into account and this figure serves as a visual reminder 
of the process. 

Finally, we mention the sensitivity of the extrapolations to the range in iV max used to fit 
the exponential function. Results at iV max = 10 with NN-only allow an assessment of the 
extrapolation procedure. We find that including the iV max = 2 points bias the extrapolation 
high, so the best estimates use iV max = 4-8 when NNN matrix elements are included. 

Figure M shows the convergence of the 6 Li ground-state as a function of Ml for selected 
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FIG. 6: (Color online) Ground-state energy of 6 Li as a function of the oscillator frequency, Ml. 



12, [32] up to N A3r 



40. 



Again we use a N 3 LO NN [29(] up to iVA2max = 300 and N 2 LO NNN 
The left and center panels are NN-only while the right panels include both induced and initial 
NNN. The black arrows indicate the experimental value. Note for comparison, the left is shifted 
relative to the center due to the difference in initial interactions. 



A's. The separate panels compare, for two different A's (top and bottom), the current results 
with a previous study [2| of NN-only calculations where the evolutions were performed in 
momentum space. The momentum-space evolutions used only the neutron-proton part of 
the interaction, V np , as an average for the complete NN interaction and added the Coulomb 
contribution separately afterward, causing a systematic overbinding of about 1 MeV. We 



have used the whole N 3 LO 29] interaction and included the Coulomb in our evolutions, but 



have checked that we recover the previous NN-only results with V np . The left panel shows 
the momentum-space evolved calculations. The center panel shows a reproduction of those 
results with a systematic shift due to the revised handling of the initial interaction. The right 
panel shows the full NN+NNN calculation. All three panels show good correspondence in 
Ml dependence between the previous and current calculations, indicating that the NN-only 
calculations are good predictors of the minima for the larger 3N-inclusive version. 

As discussed above, a lower Ml results in a lower cutoff, A uv , (see Eq. (E}) and requires a 
larger basis (larger A^ max ) to achieve the same convergence. This can be seen in the NN+NNN 
panel by observing the trend in A" max for each Ml. The values that each of the Ml spaces are 
converging towards are different, indicating that the initial Hamiltonian has been truncated 
by different incomplete bases and will not obtain the same result; the lowest N max curve is 
not yet flat. This is especially evident at the smallest Ml which has the worst truncation. 
Improving this convergence requires increasing the three-body basis size, A^max- Here, at 
^A3max = 40, the dimension of a single A = 3 channel (with quantum numbers J n T) can be 
7-8000 states, requiring 60-70 nodes for 12 hours with a hybridized MPI-OpenMP evolution 
code. This is currently the most significant computational bottleneck in our SRG program. 

In Fig. [7J we show the spectrum for 6 Li in absolute level energies on the left and excited 
state spacings on the right. We chose A = 1.8 fm -1 due to sufficient convergence as indicated 
by our extrapolation procedure. This example happens to closely match the excitation 
spectrum of available LS based results, though the variation in A is not large. We also 
include a spectrum with the excitation energies shifted to the converged ground state energy. 

The upper pair of plots shows the convergence with respect to the final A" max of the 
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FIG. 7: (Color online) Spectrum of 6 Li as a function of ^ASmax- On top (bottom) shows convergence 
in N max (A^A3max)- Left (Right) show absolute (excitation) energies. All four panels are for A = 
1.8 ml- 1 with NN+NNN at HQ = 20 MeV and iY A 2ma X = 300. 



6 Li calculation. In the upper left pane we can see a consistent convergence pattern from 
2hQ down to 8MI, but the results are clearly not converged at 8M1 On the right, the 
excitation energies indicate that the higher J states are converging more slowly due to 
stronger dependence on the higher A=3 partial-waves. This was tested by using various 
levels of truncation in the initial Hamiltonian for channels with higher values of J so that, 
for example, if J > 1 then iV max < A^A3ma X for that channel. We can also see here that the 
excitation energies compare well to the existing LS calculations. 

The lower panes show the dependence on the initial truncation of the A=3 sector, A^A3max- 
On the left we see that the Hamiltonian is well converged when iVA3ma X = 40. On the right 
the excitation energy of the 3 + state drops significantly with increasing iVA3 max - In these 
calculations, not only is A^max larger, but individual channels are truncated less severely 
as well to provide more of the Hamiltonian for these higher J states. Specifically, at A^max 
= 36 each step up in J corresponded to a reduction in A^ max of 4. But, for the AA3ma X = 40 
calculations this was changed to 2 A r max for every step in J. It is possible to push A^max 
and these channel truncations even higher and may be needed in the future for increased 
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accuracy. 



TABLE III: Results for binding energy in all calculations for 3 H and 4 He. All errors for SRG 
results are convergence margins at the quoted basis size. The columns marked "NN+3N" show 
NN+NNN-induced values except for rows that include an initial NNN. Basis sizes are A^A3max = 
40 and A^max = 196 for all except N 3 LO calculations used AA3max = 32. 



Nucleus/ 


nn 


NN-only 


NN+3N 


NN-only 


NN+3N 


LS 


Potential 




A = 1.0 


A = 1.0 


A = 1.8 


A = 1.8 




3 H - avl8 


28/52 


-7.487(1) 


-7.486(40) 


-8.467(1) 


-7.486(40) 


-7.62(0) 


cdbonn 


28/52 


-7.505(1) 


-7.863(40) 


-8.553(1) 


-7.863(40) 


-7.99(1) 


n31o 


20 


-7.471(2) 


-7.852(5) 


-8.351(1) 


-7.852(5) 


-7.852(5) 


+NNN 


20 




-8.473(5) 




-8.473(5) 


-8.473(5) 


4 He - avl8 


44 


-24.419(23) 


-24.339(14) 


-29.267(15) 


-23.904(25) 


-24.23(1) 


cdbonn 


44 


-24.484(16) 


-26.217(9) 


-29.739(11) 


-25.926(17) 


-26.1(1) 


n31o 


28 


-24.284(0) 


-25.641(1) 


-28.446(1) 


-25.325(1) 


-25.39(1) 


+NNN 


28 




-28.661(3) 




-28.464(2) 


-28.50(2) 



Table HLT1 gives a summary of the levels of convergence achieved in the present calculations. 
These are unextrapolated results from complete model spaces for the purposes of comparison 
to existing and future results and experiment. Extrapolated results for 6 Li will be given 
below. We strongly advise the reader that the specific choice of A is less important than the 
A dependence in the final results. The A dependence is the indicator that many-body forces 
are being induced to account for the renormalized components. Here we choose to display 
A's = 1.0 and 1.8 fm _1 as they reach over the range of A dependence in this work. For the 
potentials AV18 and CD-Bonn the value for NN+NNN-induced in the table corresponds to 
the unevolved minimum. The optimal frequency for the unevolved potential in these cases 
is Ml ~ 52 MeV and these values are quoted in the table. 

IV. EVOLUTION OF MANY-BODY FORCES 

In Fig. [HJ the general flow of an SRG evolved nuclear Hamiltonian is illustrated. The 
ground-state energy of the triton is plotted as a function of the flow parameter A from oo, 
which is the initial (or "bare") interaction, toward A = 0. We used A^max = 196, A^max 
= 36, and HQ = 20 MeV, for which all energies are converged to better than lOkeV. The 
previous work ^ used a less stringent A^max = ^A3max = 36. However, those results are 
within 1 keV of the current calculations, showing that the larger AA2max is not critical for 
3 H. 

We first consider the NN-only curve (squares). If H\ is evolved in only an A = 2 system, 
higher-body induced pieces are not included. The resulting energy calculations will be only 
approximately unitary for A > 2 and the ground-state energy will vary with A (squares). 
Keeping the induced NNN matrix elements, by performing the A = 3 evolution, yields a flat 
line (circles), which confirms an exactly unitary transformation in that sector. The line is 
equally flat if an initial NNN is included (diamonds). Note that the net induced three-body 
is comparable to the initial NNN contribution and thus is of natural size. 
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FIG. 8: (Color online) Ground-state energy of 3 H as a function of the SRG evolution parameter, 
A. See Table U for the nomenclature of the curves. The calculations use A^max = 196, A^max = 
36, and Ml = 20 MeV. 



The shape of the NN-only curve, first increasing binding and then rebounding, can be 
understood from the fact that early in the evolution the high-momentum matrix elements 
are most affected by the transformations. These are predominantly the short-range repulsive 
parts of the potential and are transformed into strength in the induced NNN matrix elements. 
As the evolution progresses to lower momentum scales the more attractive parts of the 
potential are affected, causing a rebound in the NN-only result. Note that all this information 
is not lost but merely reorganized into NNN terms, and when we keep those terms we 
regain the unitary result. This analysis is supported by studies in one-dimensional models 
that showed purely attractive initial potentials to have monotonically decreasing binding 
energy [l3[ 21 



In Fig. [91 we examine the SRG evolution in A for 4 He using a chiral N 3 LO potential |2i 
with HQ = 36 and 28 MeV, the dashed and solid sets respectively. The (V A ) and (V^ ) 
matrix elements were evolved with basis sizes N^max = 196 and N^ max = 32 and then 
truncated to iV max = 18 at each A to diagonalize 4 He. Again, the higher A^max has little 
impact on the final results for this nucleus. The NN-only curve has the characteristic shape 
discussed above. When the induced NNN is included, the A dependence is significantly 
reduced. The pattern only depends slightly on the inclusion of initial NNN interaction. In 
both cases the dotted line represents the converged value for the initial Hamiltonian using a 
Lee-Suzuki based procedure. The residual A dependence is due to missing induced four-body 
forces. 

At large A, the discrepancy with the dotted line is due to a lack of convergence for 
unevolved potentials at N max = 18, but at A < 3fm _1 SRG decoupling takes over and the 
discrepancy is due to short-range induced four-body forces. This transition is emphasized 
by showing the calculation at two different values of HQ. The point in A where they meet 
is an indicator of the momentum scale at which the evolving Hamiltonian is converged with 
A max = 18- All the information included in the initial Hamiltonian at N^ 3mgx has been 
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X [fm *] 

FIG. 9: (Color online) Ground-state energy of 4 He as a function of the SRG evolution parameter, A. 
Note the comparison of two values for Ml (solid and dashed) and the best LS results (dotted). Here, 
A r A2max = 196 and A^max = 32 compared to A r A2max = A^max = 28 for the prior work [3:]. The 
results here are converged to within 10 keV of those from Ref. [sj]. See Table U for the nomenclature 
of the curves. Note the comparison of two values for Ml (solid and dashed) and the best LS results 
(dotted). The thick arrow indicates the experimental value. 



transformed into the smaller basis denned by iV max . Any residual difference between values 
of Ml (invisible on this scale) indicate the level of convergence with respect to the included 
(V^) matrix elements denned by N& 3max . 

In the three-body-inclusive curves the discrepancy due to induced four-body forces is 
about 50keV net at A = 2fm~ 1 . This is small compared to the rough estimate in Ref. (33[ 
that the contribution from the long-ranged part of the N 3 LO four-nucleon force to 4 He 
binding is of order a few hundred keV. If needed, we could evolve 4-body matrix elements 
in A = 4 and will do so when nuclear structure codes can accommodate them. 

Figure [10] compares the flow of 3 H and 4 He binding energies for several initial Hamilto- 



nians, AV18 [27| . CD-Bonn 28] . and N 3 LO [29|. For 3 H we have used harmonic oscillator 
frequencies Ml = 28, 28, and 20 MeV respectively. For 4 He these optimal frequencies are 
HQ = 44, 44, and 28 MeV respectively. The general shape of all the NN-only curves is quite 
similar here with a initial dip in binding and then a turn over at A = 1.8 fm -1 . Evolution to 
low A (< 2.0 fm -1 ) of different initial Hamiltonians at A = 2 produce a very similar evolved 
form [18j. Previously this had only been observed at the level of comparing selected two- 
body matrix elements. Here we show that observables at A > 2 also exhibit this behavior. 
The NN-only points for all three initial Hamiltonians converge as A decreases past 1.8. 

Note that the values for the unevolved potentials AV18 and CD-Bonn do not approach 
the Lee-Suzuki results because these potentials require a larger UV cutoff and therefore a 
larger Ml (~ 52) or iV max for convergence. The converged values for the bare potentials are 
quoted in Table II III and would be the equivalent of the unitary line shown in Fig. [SJ 

Of course, we should not expect the NN+NNN-induced calculations to produce identical 
results at small A because they are not equivalent Hamiltonians at the A = 3 level. However 



16 



-6.5 



> 



W 



00 

e 

s 

a 




X[fm ] 



FIG. 10: (Color online) Ground-state energy of 3 H and 4 He as a function of A, starting from 



Hamiltonians based on CD-Bonn 



28] 



Argonne V 18 [27|], and N 3 LO |2g|. See Table H for the 
nomenclature of the curves. Here we used Ml = 28/28/20 and 44/44/28 for 3 H and 4 He respectively. 
The dotted lines show the converged LS results for each potential. For 3 H these values are AV18 
= -7.62, CD-Bonn = -7.99, and N 3 LO = -7.85. 



their very similar shape indicates a specific scale dependence of three- and four-body forces 
generated during evolution. This is reminiscent of evolution of the chiral interaction in Fig. 
with and without initial NNN, where the shape is similar but shifted by the initial difference. 
This is a promising indication of a universality phenomenon. A full test of this idea will 
require coding analogous three-body interactions for the other initial potentials (i.e., IL-IX 
for AV18) and evolving other initial NN interactions. 

Figure [11] shows extrapolated ground-state energies for 6 Li at different values of A. We 
used truncations A^max = 300 and A"A3max = 40 and performed the final diagonalization up 
to N max = 8. The gray bands represent the best Lee-Suzuki results available for NN-only 
and NN+NNN initial interactions, with error estimates. The analogous truncations for these 
calculations were A^max = 400 and A"A3max = 40 with 6 Li being calculated up to A" max = 
14 for NN-only. The results are 28.5 ± 0.5 MeV without and 32.5 ± 0.5 MeV with initial 
three-body forces. 

The A dependence is shown for the lower values where the result is near convergence. 
The results plotted here are obtained from the extrapolation procedure previously described. 
This procedure accounts for the variation, with A, of the minimum Ml and extrapolates in 
JV^ the converged binding energy. The error bars are dependent on the range in Ml for 
which we have results. For any given A the error bars will be larger if the optimal Ml is not 
present in the data set used for extrapolation. In fact the extrapolation tends to predict 
too low (more negative) as measured by the predictions of A^ max = 10 points in the NN-only 
curve. This feature has also been confirmed as we systematically added to the data set for 
the A's with the largest error bars; The extrapolated points rose with better values of Ml, 
flattening the curve and reducing the apparent A dependence. 

The hierarchy of induced many-body forces can be assessed in fig. [11] by comparing the 
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FIG. 11: (Color online) Extrapolated ground-state energy of 6 Li as a function of the SRG evolution 
parameter, A. Error bars are based on variation of the minimum in Ml as discussed in the text. 
The truncations used were Nj&max = 300 and iVA3max = 40. 



spread in the NN-only curve to that of the NN+NNN curves. To do so, note that the NN- 
only curve must coincide with the LS result at A = oo. The spread in A has been reduced 
from 4 MeV to < 1 MeV. The majority of induced many-body forces missing from the 
NN-only curve is due to three-body forces subsequently included in the other curves. Notice 
that the shape of the evolution curve is very similar to those of 4 He from any of the initial 
potentials used — a gentle rise to 1.8 and a slightly steeper slope down through 1.0. We 
interpret this to indicate that the majority of many-body (A > 4) forces induced are for 
A = 4 and that 5- and 6-body forces are significantly smaller. This is consistent with the 
expected hierarchical flow of induced many-body forces [l3[ • The spread here is roughly one 
MeV, compared to the 30-60 keV found in 4 He calculations. 

Our NN-only curve is almost identical in shape to previous momentum space studies, de- 
spite the difference in initial Hamiltonians used. The previous study used only the neutron- 
proton interaction for all NN interactions while we have used the full isospin-breaking po- 
tential. The results are systematically shifted up (in the previous study 6 Li was overbound 
by ~ lMeV) in relation to the NNN-inclusive curves. The error bars for the same iV max 
are roughly the same as the previous work, confirming that NN-only calculations are good 
predictors of appropriate Ml values. 

Note that the many-body forces do not explode as previously feared and that evolution 
to lower A may not be unreasonable. While previous NN-only studies showed induced 
three-nucleon forces growing uncontrolled below A = 1.5, we see that inclusion of these 
matrix elements produces a more gentle A dependence. So, evolving to lower A to improve 



convergence may be useful in future calculations. Recent results [25| suggest that, for the 
choice of SRG generator in Eq. (j2J), the many-body forces may grow with A, so monitoring 
in the rest of the p-shell (with adequate codes) is vital. Alternative choices for the SRG 
generator and sophisticated extrapolation techniques may play a central role. 
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TABLE IV: Results for extrapolated binding energy of Li at various values of A and includes error 
bars from the extrapolation. The analogous results from Lee-Suzuki calculations are —28.5 ± 0.5 
MeV for NN-only and NN+NNN-induced, and -32.5 ± 0.5 MeV for NN+NNN. The experimental 
value is —31.99 MeV. The HQ and range in A max used for each extrapolation is also quoted. 



A 


best h£l 


NN-only 


NN+NNN-induced 


NN+NNN 






(AW 4-10) 


(A max 4-8) 


(A max 4-8) 


2.5 


24 


-31.3 ± 0.5 


-28.5 ± 1.0 


-32.9 ± 1.0 


2.2 


20 a 


-31.6 ± 0.3 


-28.3 ± 0.6 


-32.4 ± 0.6 


2.0 


20 


-32.0 ± 0.2 


-28.1 ± 0.3 


-32.2 ± 0.3 


1.8 


16 


-32.8 ± 0.1 


-28.0 ± 0.2 


-32.2 ± 0.2 


1.5 


16 


-33.00 ± 0.05 


-28.1 ± 0.1 


-32.3 ± 0.1 


1.2 


16 


-31.85 ± 0.05 


-28.4 ± 0.1 


-32.75 ± 0.1 


1.0 


16 


-29.75 ± 0.02 


-28.8 ± 0.1 


-33.2 ± 0.1 



"This one point has a different optimal Ml for NN-only at Ml = 24. 
V. HIERARCHY 




FIG. 12: (Color online) Expectation value of the components of the nuclear interaction in the 3 H 
and 4 He ground states as a function of the SRG parameter A. The initial interaction was N 3 LO 
NN H and N 2 LO NNN at Ml = 20 and 28 MeV (left and right) and A A3m ax = 32 with the 4 He 



final truncation at A max = 18. The insets show expanded details for the 3- and 4-body forces. 



In order to more fully understand the SRG evolution we can trace the individual parts 
of the Hamiltonian. Figure [12] shows the 3 H and 4 He ground-state expectation values of 
individual components of the evolving Hamiltonian as a function of A. The insets show an 
increased scale for closer inspection of 3- and 4-body expectation values. Here a hierarchy of 
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induced many-body forces is evident. The magnitude of variation in A for each curve differs 
by approximately an order of magnitude. Cancellations between T re i and Vnn are reduced 
significantly over the course of the evolution. The strength of matrix elements at large 
momenta is being reorganized into (shifted to) lower momentum matrix elements. Hence 
the absolute values of the expectation values (T rel ) and (Vnn) are reduced. 

Also, note the correspondence between the (V^) curve and NN-only evolution curves 
such as in Figs. [5] and [HI The size of the three-body force reaches a minimum corresponding 
to the point (A ~ 1.8 fm^ 1 ) of maximum binding achieved by the NN-only calculations. This 
is simply the explicit plotting of the many-body forces that is implied by the approximately 
unitary curves shown in Section IIV1 In this case, the induced three-body forces effectively 
cancel out the initial three-body terms; the expectation value, (V^n), drops almost to zero. 

To make a connection between the individual terms in the three-body interaction evo- 
lution and the running of the ground-state energy, we need the evolution equations for the 

(3) 

expectation value of K in the ground state. Denoting the ground-state wave function for 
the A-particle system by \ipf), it evolves according to (it is convenient here to use the flow 
parameter s = 1/A 4 ) 

I^HWio), = , (9) 

where U s is the SRG unitary transformation at s and 

dU, 



Vs 



u\ = -4 . (io) 



ds 

Then the matrix element of an operator O s evolves according to 

^MOsHt) = (vtfl^ - h,oM) • (ii) 

If the operator O s is transformed as O s = U s O s= oUj, then the matrix element on the right- 
hand-side of Eq. ( 11 ip vanishes, as when O s = H s . 

However, if we wish to see how one part of H s evolves, such as the expectation value of 
V^ 3 \ we obtain 

^ik (3 w> = <^ a i^ - h, vm?) , (12) 

which does not give zero in general because v} 3 ^ ^ U s V}!j U%. In the two-particle case, the 
analog of Eq. ( fl2|) gives d{V^)/ds = T re i]). In the three-particle case, we can expand 
Eq. (fT2|) as 



^?\v} 3) \i>f) = (^\[Vs,H s ] 3 -[ Vs ,V s {3) Mf) 

= <^lF?U cl ] + [Fi 2) ,K (2) ] c + Fi 2) ,K (3) ] + Ff,K {2) ] + Fi 3) ,K {3) ] 

-Fi 2) ,K (3) ]-Fl 3) ,K (3) ]l^) 
= mvf\H s ] + [Vf\v}%-[vf\v s ^m 

= (^\[vf\v}\-[vf\vmf), (13) 

where vf 1 and are the commutators V a 2 = [T rcl , V^] and V"/ = [T ie i, v} 3 ^}. In the 
third line, the expectation value of the commutator, [V s ,H S ], vanishes identically. 
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The subscript "c" in the first term indicates that only connected parts of this commutator 
have been kept, and refers to a diagrammatic formalism developed in Ref. [19). Computing 

[V s , V^} in the three-particle space involves all nucleons democratically. However, com- 
mutators which leave one nucleon as a spectator cancel out in the A = 2 sector. So, we must 

compute [V s , V^] for the A = 2 sector and embed it in A = 3 so we can isolate the piece 
that affects the evolution of V^ 3 \ In general, this subtraction is required at all sectors in A, 
and the "c" here indicates that this procedure has been done. 
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FIG. 13: (Color online) Contributions from individual terms to d(V^ )/d\, the flow of the triton 
ground-state expectation value of the three-body potential. On the left (right) the calculation is 
done without (with) an initial three-body interaction. 



In Fig. H3] we show the ground-state expectation values of the terms in Eq. (1131) for the 
triton. The left panel shows the calculations with just induced NNN interactions and the 
right panel with an initial three-body force. It is most useful for our analysis to convert 
from derivatives with respect to s to derivatives with respect to A using -j- = — ^r-yr- The 
dominant contribution to the evolution of the three-body potential matrix element is the 
two-body connected part, [V s , Vs ) c - This dominance is stronger here than seen in a one 
dimensional analog [13| perhaps due to a stronger initial hierarchy in the EFT compared to 
the initial conditions chosen in one-dimension. Again, the evolution of three-body matrix 
elements depends on the interplay between long- and short-range, attractive and repulsive 
parts, which lead to scale-dependent inflection points and slopes. 

We can repeat the above analysis for A = 4 and obtain 

^(^ 4) IK (4) I^ 4) ) = (4nvf\vj% + [Vf\v s (% 

+Ff\v s (%-Ff,v^]m, (14) 

where we find no fully connected terms with only two-body forces. Again, disconnected 
terms involving two and three body potentials cancel out in the lower sectors. The leading 

(2) (3) 

terms are commutators with one V s " and one Vs , followed by connected terms quadratic 
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Ill V s w and one term quadratic in V} . All terms are small and additional cancellations 
among them may further suppress the four-body contribution. Thus, the initial hierarchy of 
many-body forces suggests that induced four-body (and higher-body) forces will be small. 




FIG. 14: (Color online) Contributions of the terms in Eq. Q3] to d(V^ }/d\, the evolving 4 He 
ground-state expectation value of the four-body force. On the left (right) the calculation is done 
without (with) an initial three-body interaction. 



In Fig. [T3] we plot these contributions to the evolution of the four-body expectation value. 
On the left panel is shown the calculations with just induced NNN interactions and the right 
panel includes initial three-body forces. Again it is more useful to convert the derivatives 
in s to derivatives in A. The interplay of contributions is much more complicated than for 

(2) (3) 

A = 3. We can see cancellations between one commutator involving V s and V s (blue 
diamonds), and the term quadratic in V s (green triangles). This is in slight contrast to the 
analogous case in one dimension where all four terms were involved in less straightforward 

(2) 

cancellations. No terms quadratic in Vs appear because no connected diagrams can be 
constructed for the four-particle evolution. The total derivative of V s ^ is small until below 
A = 2. Again the dominant contribution to the flow is the lowest order commutator, and no 
feedback in Vs^ is present. 

There is room here for dependence of the induced many-body forces on the strength of 
the initial three-body potential and figure HH supports this as far as A = 4. Forthcom- 



ing results 25] provide evidence of such dependence increasing with A. Other forms of 
SRG, such as one with the replacement T rd — >■ T rd + V 2n , may be useful in controlling the 
renormalization of the long range parts of the initial potential. 

We also note that any complete analysis of the growth of induced many-body forces must 
involve converged or extrapolated results at the optimal Ml for each A. The analysis tool 
shown here is only meaningful when viewed at a single Ml over the course of evolution in A 
and direct comparison to plots of the type shown in figure [TT1 is difficult. 
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VI. OBSERVABLES 



While accurate reproduction of nuclear binding energies is the first step in nuclear struc- 
ture calculations other observables can offer additional information about the effects of 
renormalizing high-energy degrees of freedom, short-range correlations, and other details of 
a properly fit initial Hamiltonian. While we know that the harmonic oscillator basis is not 
an ideal environment for certain long-ranged observables, such as the rms radius, we have 
existing Lee-Suzuki renormalized benchmarks with which to compare. And electromagnetic 
transitions, such as B(E2)'s and B(Ml)'s are notoriously difficult, both to calculate and to 
measure, making this an important area of prediction for theory. All such observables are an 
important next test in understanding the quality of the many-body wavefunctions resulting 
from SRG-evolved interactions. Here we present a small sampling of results, focusing on 
convergence patterns. 

Here we are plotting the unevolved operator expectation value in the evolved wavefunc- 
tion. This is a reasonable way to visualize the effect of evolution on the structure with 
respect to particular operator. However, consistent renormalization of the operators them- 
selves is an important part of a robust nuclear structure program. Work along these lines is 
proceeding and is partly presented in Ref. [3]. Extending beyond A = 2 will be covered in 
a forthcoming paper. 
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FIG. 15: (Color online) Various observables in 6 Li as a function of N, 
results are with included initial NNN forces and iVA2max = 300 and Na3 
MeV. The arrow shows the experimental value and the dotted line shows the best LS result. 



luax for select A's. These 
1X = 40 and HQ. = 20 



Shown in Fig. [15] are selected observables for 6 Li as a function of A max up to A max = 
8. Included are simple extrapolation curves shown by the dotted lines extending from the 
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data points. Table |V] shows the values for iV max = 8 at select values of A. In all cases 
the extrapolated values compare well to the established Lee-Suzuki based results, but show 
room for improvement with respect to the experimental values. Note the small scales on 
some of the plots, especially the quadrupole moment and B(M1;2 + 1 — > l + 0). 

Some of these observables exhibit a non-variational pattern in N max , such as the 
quadrupole moment and the B(M1) shown here. These operators have strong coupling 
between shells of iV max and iV max +2 that result in complex cancellations from one trunca- 
tion to another. However, SRG evolution seems to improve the variational properties of 
these observables. Access to larger iV max model spaces will facilitate further study of these 
quantities. 



TABLE V: Results for selected properties of Li. Here we choose A = 1.0 and 1.2 fm~ due to 
their convergence properties. These results were obtained in a basis space with iV max = 8 and hQ 
= 20. All results have included initial 3N forces at N 2 LO. The LS results were obtained at N max 
= 8 with Ml = 13 MeV 36]. Note that Kl values for LS and SRG procedures do not necessarily 



correspond to one another. 



Observable 


A = 1.0 


A = 1.2 


Expt. 


LS 


r p [fm] 
Q(l+0) [e fm 2 ] 
B(E2;3+0 -> 1+0) 
B(M1;0+1 1+0) 
B(M1;2+1 ->■ 1+0) 


2.2841 
-0.0132 

4.0663 
16.1499 

0.0622 


2.1913 
-0.0199 

3.8087 
15.8706 

0.0784 


2.32(3) 
-0.082(2) 
10.69(84) 
15.43(32) 
0.149(27) 


2.110 
-0.085 
3.5725 
15.0717 
0.0936 



VII. CONCLUSIONS 

We have presented ab initio calculations of several light nuclei using SRG-evolved three- 
nucleon forces. The results have smooth convergence qualities with respect to basis size, 
which enable reliable extrapolations. The extrapolated (and converged where available) 
values are within the error bars of the best existing Lee-Suzuki based calculations. Inves- 
tigating the A dependence of induced many-body forces, we find that they do not grow 
substantially as A is lowered and the range of these effects is within the established LS error 
bars. Analyzing the mechanism of flow for many-body terms reveals that the SRG is driven 
by the natural hierarchy of the initial Hamiltonian and that it preserves this hierarchy dur- 
ing evolution. This is qualitatively consistent with studies of the same in one- dimension. 
Finally we present some first results of various observables using SRG evolved many-body 
wavef unctions. 

Our results here have focused mainly on 6 Li observables and analysis in the A = 3 and 
4 sectors. However, the input Hamiltonian files produced for this work are universally valid 
for further calculations in the p-shell nuclei. Here, we were limited in basis size (to N max 
= 8 in 6 Li), but plan to apply the evolved potentials at larger A using codes capable of 
larger basis sizes. We are first interested in studies of 8 Be, 10 B, and 12 C, but this list will 
undoubtedly expand. Also, we hope these potentials will be applied using coupled cluster 



methods for even larger A [38J , and look forward to applications of SRG evolution to external 
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operators. Our work here provides no indications of problems as high as 6 Li with T re \ as the 
SRG generator. Other forms of the SRG generator may be useful in controlling the growth 
of many-body forces in other nuclei. 

In addition to the above ongoing work, we will apply the evolved three-body interactions 
developed here to NCSM/RGM calculations [39] of light nuclear reactions. The NN-only 
evolved interactions have so far produced good scattering and reaction results for s- and light 
p-shell nuclei. Adding the evolved three-body interaction to the NCSM/RGM formalism will 
further improve accuracy and allow us to extend its applicability to heavier p-shell and light 
sd-shell nuclei. 
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